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model taking into account effects of spatial non- locality 

PACS codes: 02.30.Jr; 47.50. Cd; 83.10. Gr 



(N 

Oh 

_^ [ Keywords: Wave patterns; compactons hydrodynamic-type model; spatial and 

temporal non-localities 

"> ■ 1 Introduction 

(N 

w I In this paper evolutionary PDEs describing wave patterns with compact support are 

fSj ■ studied. Very often wave patterns play key rules in nonlinear transport phenomena 

t:j- ■ [IIISIE]. One of the most advanced mathematical theory dealing with wave patterns' 

O , formation and evolution is the soliton theory. The well-known Korteveg - de Vries 

(^ [ equation (KdV) 

Ut + l3uUx + Uxxx = 0, (1) 

S^ . possesses a one-parameter family of solutions describing exponentially localized wave 

12 a^ 
u{t, x) = ——sech'^[a{x -Aa^t)]. (2) 



Solitons demonstrate many outstanding features. Being the solutions to non-linear 
evolutionary equation (Q, they manifest "elastic properties" during the collisions. 
Besides, any sufficiently smooth initial data possessing a finite energy norm creates 
a chain of solitons moving with different velocities. All this is recognized as a 
consequence of the complete integrability of equation ([T]) and existence of infinite 
hierarchy of conservation laws. 

Actually the solitons have been associated with a number of phenomena observed 
in science and technics. However, solution ([2]) cannot be fully identified with any 
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kind of solitary waves occurring in natural phenomena since the wave pattern de- 
scribed by this formula at fixed t extends to ±00. 

In recent years there have been discovered another type of solitary waves sup- 
ported by the following generalization of the KdV hierarchy [5]: 

Kim, n) = u, + (un, + («")... = 0, m, n > 2. (3) 

In case when m = 2, n = 2 a. generalized solution to this equation is described by 
the following formula |5J: 

{+ ^ _ / ^ ^0^2 [{x -Dt)/A] when \x-Dt\<27i, , . 

^^'^^~\ when |a;-Dt| > 27r. ^^ 

Solutions of this sort are rather typical to the whole K{m, n) hierarchy. All of them 
vanish outside some compact domain and hence they are referred to as compactons. 
Study of K{m, n) hierarchy is actually in progress and there has already been real- 
ized that most of these equations are not completely integrable and do not possess 
an infinite set of conservation laws. Nevertheless, compacton-supporting equations 
inherit many features of equations belonging to the KdV hierarchy. In particular, a 
sufficiently smooth perturbations taken as Cauchy data give rise to a chain of com- 
pactons O El [7l [8] . Besides, compactons manifest almost perfectly elastic features 
during the mutual collisions [HI [7] . 

It should be emphasized, yet, that most papers dealing with the subject in ques- 
tion, are concerned with the compactons being the solutions to either completely 
integrable equations, or those which produce a completely integrable ones when 
being reduced onto subset of a traveling wave (TW) solutions [5], [9], |10] . 

In this paper compacton-like solutions to the hydrodynamic-type model taking 
into account the effects of spatial non-locality are considered. The model occurs to 
share the compacton-like solutions with a Hamiltonian system but only accidentally, 
i.e. at certain values of the parameters. In spite of such restriction, existence of this 
type of solutions seems to be significant for their appearance is connected with the 
presence of non-local effects and rather cannot be manifested in any local hydrody- 
namic model. Besides, localized invariant solutions sometimes manifest attractive 
features and when this is the case, they can be treated as some universal mechanism 
of the energy transfer in media with internal structure leading to the given type of 
the hydrodynamic-type system. 

The structure of the paper is following. In section 2 we introduce dynamic equa- 
tions of state (DES) taking into account effects of spatial and temporal non-localities. 
In section 3 we perform the qualitative analysis of a hydrodynamic system of balance 
equations closed by DES taking into account effects of spatial non-locality and state 
the conditions leading to the existence of compacton-like TW solutions. In the last 
section we summarize the results obtained, discuss open problems and prospects of 
further developing the actual study. 



2 Non-local hydrodynamic-type models 

We are going to analyze the existence of compacton-like solutions within the hy- 
drodynamic-type models taking into account non-local effects. These effects are 
manifested when an intens pulse loading (impact, explosion, etc.) is applied to media 
possessing internal structure on mesoscale. Description of the non-linear waves 
propagation in such media depends in essential way on the ration of a characteristic 
size d of elements of the medium structure to the characteristic length A of the wave 
pack. If d/\ is of the order of 0(1) then the basic concepts of continuum mechanics 
are not applicable and one should use the description based e.g. on the element 
dynamics methods. The models studied in this paper apply when the ratio d/\ is 
much less than unity and therefore the continual approach is still valid, but it is not 
as small that we can ignore the presence of internal structure. 

As it is shown in [H], in the long wave approximation the balance equations for 
momentum and mass retain their classical form, which in the one-dimensional case 
can be written as follows: 

«t+Px = 0, (5) 

Pt + p'Mx. = 0. (6) 

Here u is mass velocity, p is density, p is pressure, t is time, x is mass (Lagrangean) 
coordinate. So the whole information about the presence of structure in this ap- 
proximation is contained in DES which should be added to system (l5])-(l6]) in order 
to make it closed. 

It is of common knowledge that a unified dynamic equation of state well enough 
describing the behavior of condensed media in a wide range of changes of pressure, 
density and regimes of load (unload) actually does not exist. Various particular 
equations of state describing dynamical behavior of structural media are derived by 
means of different techniques. There exist, for example, several generally accepted 
DES in mechanics of heterogeneous media, derived on pure mechanical ground (see 
e.g. IHIITS]). Contrary, in papers [HI HH [16], [1^ derivations of DES are based on the 
non-equilibrium thermodynamics methods. Both, mechanical and thermodynamical 
approaches under the resembling assumptions give rise to similar DES, stating the 
functional dependencies between p, p and their partial derivatives. 

There is also a number of works aimed at deriving the DES on the basis of the 
statistical theory of irreversible processes (see [HI [19] and references therein). It 
is rather firmly stated within this approach that DES for complicated condensed 
matter, being far from the state of thermodynamic equilibrium, takes the form 
of integral equations, linking together generalized thermodynamical fluxes /„ and 
generalized thermodynamical forces L^, causing these fluxes: 

/t r+oo 

dt' K^n{t,t';X, X') Qra [Lk {f , x')] . (7) 

'OO J —oo 

Here -ft'mn (^? i'; ^i x') is the kernel of non-locality, which can be calculated, in prin- 
ciple, by solving dynamic problem of structure's elements interaction. Yet, such 



calculations are extremely difficult and very seldom are seen through to the end. 
Therefore we follow a common practice [201 HI] cind use some model kernels describ- 
ing well enough the main properties of the non-local effects, in particular, the fact 
that these effects vanish rapidly as |t — t'| and |a; — a;'| grow. 

DES derived here are based on the following relation between the pressure and 
density: 

p{t, x) = f [p{t, x)]+ \ K{t-t', x~x') g[p{t', x')] dx'\ dt'. (8) 



Let us ffist assume that effects of spatial non- locality are unimportant. In this case 
the kernel of non-locality can be presented as i^ (t — t') 6 {x — x') and the flux-force 
relation (El) takes the following form: 

p{t, x) = f[p {t, x)]+ f K{t- t') g [p {f, x)] dt'. (9) 

^— oo 

Since we rather do not want to deal with the system of integro-differential equations, 
our next step is to extract some acceptable kernels enabling to pass from (^ to 
pure differential relations. Differentiating equation Q with respect to the temporal 
variable we obtain: 

Pt = f [p] pt + ir(0) g[p\+ I Kt (t - t') g [p it', x)] dt'. (10) 



Equationl TTOl) is equivalent to a pure differential one provided that function K{z) 
satisfies equation K{z) = cK{z). In this case we can obtain from ([9D, (ITOl) the 
following differential equation: 



r{pt~f[p]pt}=rAg[p] + f[p]-p. (11) 



This equation corresponds to fading memory kernel K{z) = A exp [— -]. For A = 
1, f [p] = XP^^^i 9 [p] = — o"p"+-^, we get the following DES: 

r{pt-x{n + l)p"pa = «:p"+^ -p. (12) 

where k = x — <^t- Equations coinciding with (lT2l) under certain additional assump- 
tions are widely used to describe nonlinear processes in multi-component media with 
one relaxing process in the elements of structure fTH [TB] . 

In more complicated cases, e.g. when more than one relaxing process should be 
taken into account, another kernels are to be applied. Let us look for the kernels 
leading from Q to pure differential equation under the assumption that such passage 
is possible with the help of two differentiations. Taking derivative of flTOl) with 
respect to the temporal variable we get: 



/ [p] Pt] + i^(0) g [p] Pt + k[0] g[p]+ f Ktt {t - t') g [p (f, x)] dt'. (13) 
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So now we are interested in whether the hnear combination hpu + spt + p can be 
expressed as a pure differential equation. It is easy to see that the integral term 
corresponding to this combination is as follows: 

f {hKu {t - t') +sKt{t-t')+K{t- t')} g [p{t', x)] dt'. 



Our requirement will be addressed if K[z] satisfies the equation 

h K,, {z) + sK,{z) + K {z) = 0. 
If this is so, then DES takes on the form 

hptt + spt+p = f{p) + [sK{0) + hKt{0)] L{p) + 



(14) 



(15) 



+ 



sf{p) + hK{0)g{p) pt + h f(p)p,, + f\p)p 



Their are two cases for which requirement of fading memory holds true: 

z_ z_ 

K[z\ = a e ^1 + /3e ^2 , h = T1T2, s = n + T2, 
and 



K[z] 



e ^1 sm 



r2 



+ 7 



h 



rl + rr 






In the first case we have the kernel describing media with two relaxing processes, 
while in the second case the kernel describes relaxing media with oscillating compo- 
nent in the elements of structure. Note that in both cases we get formally identical 
DES flTSj) . yet the parameters corresponding to them satisfy quite different inequal- 
ities. In the first case the parameters /i, s obey the inequalities s^ > 4/i > 0, while 
in the second case the inequalities are as follows: < s^ < 4 /;,. So the parameters 
lie in the distinct sets of parameter space and this fact occurs to be crucial since the 
behavior of the system of balance equation ([5])- ([6]) closed by DES (fTSl) is extremely 
sensible on the values of these parameters. pLTj i22j . 

Now let us address the case of pure spatial non- locality. Following [21] , we shall 

use the kernel of the form K\t — t' , x — x'] = a exp — (^7^) ■ 5[t — t'] giving the 
equation 



/n2 



g[p{t, x')] dx. 



(16) 



P = / (p) + '5" 

Since the function exp — (^7^) extremely quickly approaches zero as |x — x'| 

grows, we use the following approximation for function g [p {t, x')] inside the inner 
integral: 



g[p{t, x')] =g[p{t, x)] + {g[p{t, x)]}^ 



^'^' ''^+{9[pit,x)]}J^^^+0{\x'-xf) 

(17) 



1! 



Dropping out the term 0{\x' — x\^) and integrating over dx', we get 

p = f[p{t,x)]+aog[p{t,x)] + a2 {g[p{t,x)]}^^ {11 



where 



i ax . ao = a 



ao = a e \- ' i dx', a2 = a j e I- ' J dx'. 

J —oo J —oo ^ • 

It is obvious, that one can easily combine two types of non-locahties by taking 
the kernel in the form of the product 



K[t-t', X- x'] = K[t-t']- exp 



, , 2 
X — X 



I 



with K [t — t'] satisfying JHM . Employment of technics used in the last two items 
enables to obtain the DES taking into account both effects of spatial and temporal 
non-localities: 

hptt + spt + p = f[p] + [sK[0] + hKt[0]] [aog[p] + a^ {g[p] p.) J + (19) 
+s f[p] pt + h K[0] [ao g[p] + as ig[p] p.) Jt + h [f[p]pt] ^ ■ 

Let us note that DES with higher derivatives were employed in many papers 
dealing with the structured media. Qualitative and numerical investigations un- 
dertaken in [221 ESI 121] show that the modeling system taking into account both 
types of non-localities possesses a rich set of the TW solutions, including periodic, 
quasi-periodic, multi-periodic, soliton-like and chaotic regimes. 

3 Compacton-like solutions within the hydrody- 
namic-type model taking into account spatial 
non-localitiy 

Before we start to "capture" compacton-like solutions within the hydrodynamic 
model, it is desired to get a clear idea on what is a compactons from geometric 
point of view. In paper ^25] such analysis has been performed for the equations 
belonging to Rosenau-Hyman K{n, m) hierarchy and this is the way things stand. 
Being treated as function of a single value C, — x — Dt, solution (jlj) satisfies some 
second-order ODE which is obtained from ([3]) by means of the ansatz u{t, x) = U{^). 
Substituting this ansatz into the source equation one gets, after one integration and 
some manipulation, a second-order Hamiltonian system [25]. This system occurs 
to have a one-parameter family of periodic phase trajectories and the homoclinic 
trajectory surrounding them. The last one just corresponds to compacton. 

In the case of KdV the similar procedure leads to the Hamiltoinan system with 
the same geometry of the phase trajectories. The difference between two phase 



portraits arises from the fact that in the latter case the honiochnic loop is the 
bi-asymptotic trajectory of a simple saddle while the trajectory representing the 
compacton is the separatrix of a saddle settled on the line of singular points of the 
vector field. Therefor the homoclinic trajectlry is penetrated in a finite "time" and 
the compacton solution (jll) consists of non-zero part with compact support, glued 
with zero solution corresponding to the saddle point. 

Let us note at the end of this short introduction, that we do not distinguish 
solutions having a compact support and those which can be made so by proper 
change of variables. 

Now let us consider the system ([5]), ([6]) closed by DES 

p = /[p] + «:p"+2 + a [p^+Vx], (20) 

corresponding to pure spatial non-locality. Here p = p — po, where < po is a 
constant, f[p] is a function which will be defined later on. We use the ansatz 

u{t, x) = U{x-Dt) = f/(0, p(t, x) = R{0, (21) 

enabling to factorize system ([6]), ([5]), closed by the DES (!20l) . Inserting (121]) into 
equation (Q, we get the following quadrature: 

Po R[Q + Po 

Constant of integration have been chosen in such a way that u{t, ± oo) = 0. 

Using the ansatz (pTI) . we express DES in new invariant variables. Inserting ( 120|) 
into ([S]), using the formula (B^ and integrating once the expression obtained this 
way we pass, after some manipulation, to the following second order ODE 

-^— + f[R] + kR^+^ + a [i?'"+i R'Y = E = — + f[0]. (23) 

R + Po Po 

It is obvious that above equation can be re-written as an equivalent dynamical 
system. To do this, we define a new function W = —R'. Next we introduce new 
independent variable T such that jf = a R^^^ ip[R\j^, where (/^[i?] is an integrating 
factor which is incorporated in order to make the system Hamiltonian. With this 
notation we get the following dynamical system equivalent to 




—aifi 

ip[R] 



R\Rn+lW = -^^^^ 



D^R 



dH[R,W] 
dR ■ 



a{n + l)R- W^ + f{R) - /(O) + kR^^' - ^^^^^^^^ 

(24) 
Solving the first equation of system fl2^ with respect to function H[R, W] we obtain: 

H[R, W] = aip[R] i?"+'^ + 9[R]. (25) 



Now, comparing the RHS of second equation of system (^^ with partial derivative 
of fl25|) with respect to R, we get the system 

Rip'[R] = {n + l)ip[R], 
e'[R] = y,[R] IfiR) - /(O) + «:i?"+2 D'R 



po {R + Po) j ' 



The first equation is satisfied by the function (f[R] = C R^^^ . For convenience we 
put C = 2. Inserting (p[R] in the second equation we obtain the quadrature 

0(R) = 2/ R"» [m - /(O) + .fl- - ^-gf-^] ,R. 

So the Hamiltonian function is finally expressed as follows: 

H[R, W] = a/?2(«+i)H/2 ^ 2 /" i?"+i !f(R) - /(O) + fi:i?"+2 ^-^--\ dR. 

J { Po (-K + Po) J 

(26) 

The Hamiltonian function fl26|) is a first integral of the system fl24l) . so every phase 
trajectory can be presented in the form 



W"^ = ^^ iK -2 I i?"+^ 



f{R) - /(O) + kR^+' - ^'^ 



dR 



po {R + Po) 

(27) 
where iT is the constant value of the Hamiltonian on a particular trajectory. Now 
let us analyze formula fl27|) . If we want to have a closed trajectory approaching the 
origin, then we must properly choose the constant K and "suppress" all the singular 
terms by the proper choice of function f{R)- It can be easily shown by induction 
that the following decomposition for the last term inside the integral takes place: 

pn+2 n+2 

-/?"+^-po^"+-+(-i)''Po^"'''"''+-+(-i)"+vr'+(-i)"+' 



R + Po R + Po 

Hence 



^ |k-2| R-^' [fiR)-m + KR--^'] dR- 

2 / 73n+2 R"+l \ ^ 

- -—7^-Po^ + ■ + {-PoT^'R+i-por-^'logiR + po) . 
\n + 2 n + 1 J ) 



aR^( 
2D'' f R 



Po 

From this we conclude that the last term in (!27j) always produces singularities. More 
precisely, singularities are connected with monomials -R™ when m < 2{n + 1) and 
with the logarithmic term. Therefore the last term in 0271) should be rather removed 
by the proper choice of /(-R). 

A simple analysis shows that function 

f{R) = fo+ .t^ , + 91 R""^' + 92 i?"+' 
Po [R + po) 
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with A > will suppress singularity provided that D = ±v A. In fact, in this case 



iy2 _ I J J^ _ ^2(n+l) 



2giR , hR 



+ 



2n + 3 n + 2 



(28) 



where g2 = 92 + i^- With such a choice the only trajectory approaching the origin 
corresponds to -ft' = 0. If in addition gi = —ai < and (72 = "2 > 0, i.e. 

fiR) = fo+ .t^ , - «i i?"+' + «2 i?"+' (29) 

po [R + Po) 

then 



W = ±^J^^R-^^R- (30) 

and there is the point R^ = J" g|°^ in which the trajectory intersects the horizontal 
axis. 

In fact, under the above assumptions we get the geometry which is similar to that 
obtained when the member of K{n, m) hierarchy is reduced to an ODE describing 
the set of TW solutions [25J. To show that, let us consider the system arising from 
(!23|) under the above assumptions: 

^ = -2aiyi?2("+i) (31) 

^ = 2{n + l)aW^ i?2-+i + 2 «2 /?'"+' - 2 «i R'^^^^. 

System ( 13T1) possesses two stationary points lying on the horizontal axis: the point 
(0, 0), and the point (i?i, 0), where Ri = ai/a2. Analysis of the Jacobi matrix 
shows that (i?i, 0) is a center. Moreover, Ri < R^, when n > —2 so the closed 
trajectory corresponding to K = encircles the domain filled with the periodic 
trajectories. Using the asymptotic decomposition of the solution represented by the 
homoclinic trajectory near the origin, it can be shown that the trajectory reaches 
the stationary point (0, 0) in finite "time" so it does correspond to the compacton. 
Besides, we can get more direct evidence of the existence of compacton-like solu- 
tions by integrating equation fl30|) . In fact, taking in mind that W = —dR/d^, we 
obtain the following equation: 

1 dR 

7 



where 7 = ^^^, ^ = \h^^- Solving this equation we get: 

R\e] = l 7[l + sin(5(e-^o))] when-7r< 2,5(^-^o)< Svr 
*■ ■' 1^ otherwise. 

Using the formulae fl2(l - fl2^ and (152]) . one can easily reproduce the functions 
u, p and p. 



4 Discussion 

In this work we have shown that hydrodynaniic system of balance equations ([5])- ([6]) 
closed by DES ( l20|) possesses the compacton-like solutions. In contrast to analogous 
solutions to most of the compacton-supporting equations, the presented solutions do 
not form a one-parameter family. More precisely, for fixed values of the parameters 
A, «! and 02 in (129|) there exists exactly one pair of compactons moving with 
velocity vA in the opposite directions. Whether these solutions are of interest from 
the point of view of applications or not, depends on their stability and behavior 
during the mutual collisions. Discussion of these topics goes beyond the scope of 
this paper. Let us mention, however, that the mere fact of existence of compactons 
is the consequence of the non-local effects incorporation. To our best knowledge, 
this type of solutions does not exist in any local hydrodynamic-type model, i.e. the 
system of balance equations (jSJ-dHD closed by the functional state equation p = 
$ (p) . Let us note that invariant TW solutions very often play rule of intermediate 
asymptotics [261 122], attracting near-by, not necessarily invariant, solutions. This 
feature demonstrate compacton-like solutions of another non-local model obtained 
when the system of balance equations ([5])- ([6]) is closed by the DES (TT2|) . describing 
relaxing media [23]. Solutions with compact supports appear in this model merely 
in presence of the mass force. Exactly one compacton-like solution occurs to exist 
for the given set of the parameters, yet this solution serves as an attractor for the 
wave packs created by the wide class of the initial value problems [25] . In contrast 
to the above mentioned relaxing model, incorporation of the effects of spatial non- 
locality leads to the existence of a pair of compacton-like solutions in absence of an 
external force. In fact, it is shown for the first time that a non-local hydrodynamic- 
type model possesses more than one compacton-like solution and from now on there 
exists the opportunity to investigate their behavior during the collisions. 

There are some evidences in favor of the stable behavior of these compactons 
and their elastic collisions. The first is connected with the fact that by proper 
choice of parameters A, ai and 0:2 the evolutionarity conditions [2B| dp/ dp > 0, 
d^ p/d^p > can be fulfilled for DES fl20|) . Besides, the system of balance equations 
(IS])-© closed by the DES fl20l) possesses at least two conservation laws. Though 
actually it is not known how many conserved quantities assure the stability of the 
wave patterns during the collisions, it is almost certain that this number should 
not be infinite. For example, the members of the Rosenau-Hyman K{n, m) hier- 
archy, possessing with certain four conserved quantities, collide almost elastically 
and so is with the compacton-compacton and kovaton-compacton interactions in 
the Pikovsky-Rosenau model [7]. 

Another open question connected with the non-local hydrodynamic-type models 
is as follows. It is unknown as yet whether the compacton-like solutions exist within 
the hydrodynamic-type models closed by the DES fll9p taking into account both 
the effects of spatial and temporal non-localities. Our preliminary analysis shows 
that incorporation of terms connected with temporal non-locality makes the model 
dissipative. So the Hamiltonian formalism is of no use in this case and investigations 
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of the compacton-like regimes appearance should be based on some other technics. 
Study of these questions is actually in progress. 
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